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_^ ' Abstract. The cross sections of the pd — ^■"'HeT;, pd -^^HeTr^Tfirfi and pd -^^HeTT^TT^n^^ reactions have been mea- 

^ ' sured at the beam kinetic energies Tp = 1360 MeV and Tp = 1450 MeV using the CELSIUSAVASA detector setup. 

' At both energies, the differential cross section ^ of the 77 meson in the pd — ^-^Hezy reaction shows a strong forward- 

backward asymmetry in the CMS. The ratio between the pd -^^Hen^n^n^^ and pd ^-^HeTr'W" cross sections has 
been analysed in terms of isospin amplitudes. The reconstructed invariant mass distributions of the nn, ^Hen and ■'He2;r 
systems provide hints on the role of nucleon resonances in the 3n production process. 



PACS. 13.75-n Hadron-induced low- and intermediate energy reactions and scattering - 14.40-Be Light mesons 
25.40 Ve Nuclear reactions above meson production threshold 
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C3 1 Introduction to describe the data. The angular distributions from Refs. IHIHI 

have a strong forward-backward asymmetry with a backward 

The pd^'Hc + X reaction has long been used to study the suPPression, a maximum at cos 9^ « 0.5 and a forward plateau 

production of charged and neutral mesons and mesonic sys- f '^'P;^,^' Jf htly overlapping excess energies, there are data 

tems. Studying reactions with ^He in the final state gives in- c". Tfx l^^' ? f 

sight in the reaction dynamics involving three nucleons and in PROMICE/WASA and ANKE results. At high energies (Q > 

meson-nucleon final state interactions. 120 MeV), the data bank is scarce. Backward production of 

™ , -11 r ■ T ■ V mesons in pd — > He;? was studied at 17 different beam en- 

The pd — ?> He;? reaction has been or particular interest. ' . , V,r.„„ , rm 1^ ^ 1 

e 1 ; J- 1 ■ u UM o o .1 CI u ergies at the SPES IV spectrometer Dill . Parts or the 77 angu- 

Several studies near the kinematic threshold 1 ,2 3 ,4,5 , where „ — c^L^rT-,-,-, 

, • , . • , J • , lar distribution at T,, = 1450 MeV was measured by SPES III 

mostly .s-waves are involved in the production process, show a ^^-^ „, ^r-i oTiromrA o a hi. u \^ . j- j 

, / , rr> • > > , ■ . .J 0121. The CELSIUSAVASA collaboration has recently studied 



threshold enhancement. This enhancement has been interpreted 

. J. ^. „ . , J ^1,, , . ^ r^i »j Lne pa — ?• ne;? reacuon ai lwo oeam energies, i.e. 1 „ 

as an indication of a quasi-bound mtr^ nuclear state |6 1. Mea- ^^J, ^ ^^^^.^^^ .„....„„„^ L 

surements of the angular distribution at slightly higher en 

ergies from PROMICE/WASA [U and ANKE |8| indicate the 



the pd — ?• ^He;; reaction at two beam energies, i.e. rp=1450 
MeV and 7^=1360 MeV, which correspond to excess energies 
of 252 MeV and 299 MeV, respectively. The differential cross 



„ ^ UAnj u i section was measured in the backward hemisphere and at for- 

presence or p-waves at an excess energy or Q w20 MeV, while , , AT-i^cr.»jiriii j • 1 

,irv»ji7 I,- 1. .1 • J- J ward angles. At r„=1450 MeV, the backward points overlap 

at g ~ 40 MeV even higher partial waves are required in order o r ^ r r 
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with those from Ref. lfT2l . The angular distribution at 7^=1360 
MeV obtained with CELSIUS/WASA is the first measured at 
this energy. 

The direct production of three pions, i.e. pions which do 
not originate from e.g. w or 77 decay, has so far received Httle 
theoretical and experimental attention. In the isobar model dis- 
cussed in Ref. 1 13|, three-pion production should proceed via 
an excitation of one or two baryon resonances, like /I (1232) 
or the Roper A^*(1440), followed by their subsequent decays. 
Three-pion production in proton-proton collisions was studied 
at high energies 1 14,15,16] and at lower energies by CEL- 
SIUS/WASA 1 17 1. In the latter work, the ratio between o-{pp ^ 
ppn^n^n^) and a-{pp ppn^^jiPn^^) was measured and dis- 
cussed in terms of isospin amplitudes. The ratio was measured 
to be 6.3 ± 0.6 ± 1 .0 which suggests that the A^* ( 1440) An 
being the leading part of the reaction mechanism, in line with 
the isobar model presented in Ref. lfT3l . 

In the pd ^^He;7r;7r;7r case, it is straightforward to show that 
the total cross sections expressed in isospin amplitudes Mj-j^ 
are 

2 1 

cr(pd^^V{tn+7i-7i^)°^ — iMi P + -iMoP + cross terms (1) 
15 6 

o-{pd -> ^HeTT^jrV) ^ i. |2 (2) 

where Tt,„ denotes the isospin of the three pions. In the simple 
statistical approach as outlined by Fermi 1 18 1, all amplitudes in 
eq.[T]and eq.|2]are put equal and the cross terms are neglected. 
Though not justified, this simplification enables a rough com- 
parison between two channels for which no other, more realis- 
tic, model exists. The cross section ratio then becomes 

CT{pd^^mn+n-n^) _ 

If Mo is put to 0, the ratio becomes 4. In this work, the 
ratio has been estimated experimentally at 1360 MeV, which 
corresponds to an excess energy of e = 395 MeV for 37r^ and Q 
= 386 MeV for n+n^n^, and at 1450 MeV, which corresponds 
to 2 = 441 MeV for 3?:^ and Q = 432 MeV for n+n-n^. 

Multipion production is also interesting since it constitutes 
the most important background to other meson production re- 
actions like pd ^-^He?;, pd — >^He<jj, pd — >^He(?' and 
pd -^^Herin'^. 

This paper is organised as follows: in the next section, the 
reader is introduced to the CELSIUSAVASA experiment. In 
section [3] the measurement of the pd — ^^He?/ reaction is pre- 
sented and in section H] the pd ^'^He;7r*';7r*'7r*' and 
pd — S'^He;r+;r^:7!-'^ reactions are studied and compared. Finally 
the results are summarised and discussed in section|5] 



2 The CELSIUS/WASA experiment 

The measurements were carried out at the The Svedberg Labo- 
ratory in Uppsala, Sweden. The WASA detector 1 19 1 was, until 
June 2005, an integrated part of the CELSIUS storage ring. In 
the measurements presented here, a target of deuterium pellets 



II20II21II was used, designed for a 4;r detector geometry and high 
luminosity. 

The ^He ions were detected in the Forward Detector (FD) 
II22I . covering polar angles from 3° to 18°. The FD consists 
of the Window Counter (FWC) for triggering, the Proportional 
Chamber for precise angular information (FPC), the Trigger 
Hodoscope (FTH) for triggering and offline particle identifi- 
cation and the Range Hodoscope (FRH) for energy measure- 
ments, particle identification and triggering. Mesons and their 
decay products are mainly detected in the central detector (CD), 
which consists of the Plastic Scintillating Barrel (PSB), the 
Mini Drift Chamber (MDC) and the Scintillating Electromag- 
netic Calorimeter (SEC). Charged particles, mainly pions, are 
discriminated from neutral ones by their signals in the PSB, 
that also provides azimuthal angular information and covers a 
polar angular range from 24° to 159°. The momenta of charged 
particles are extracted by tracking in a magnetic field in the 
MDC. The SEC measures angles and energies of photons from 
meson decays and covers polar angles from 20° to 169°. 

A special trigger was developed to select events with ^He in 
the final state, based on the condition that ^He events give high 
energy deposit in the FWC and that hits detected by the FWC 
and the consecutive detectors FTH and FRH should match in 
the azimuthal angle. It was carefully checked in the offline 
analysis that the energy deposit thresholds were set sufficiently 
low to accept ^He ions in the full energy range, i.e. giving an 
unbiased ^He sample. 

In the offline analysis, the ^He ions are identified in the FD 
by first obtaining a preliminary particle identity (PID) using 
the AE-E-method. In short, we compare the light output in the 
detector layer where the particle stops to the light output in 
the preceding layer. The of the PID hypothesis was then 
calculated by comparing the measured energy deposits in all 
detector layers traversed by the particle to the calculated energy 
deposits. Particle hypotheses giving a larger than a certain 
maximum value were rejected. For details, see Ref. II23II24II . 



3 The pd -^^Herj reaction 

The WASA data coflected at Tp ^ 1360 MeV and Tp ^ 1450 
MeV correspond to excess energies Q = 252 MeV and Q = 
299 MeV and to 77 CM momenta of p* = 516 MeV/c and p* = 
568 MeV/c. Here and in the following, the star indicates that a 
kinematic variable is in the CM system. 

The WASA Forward Detector does not cover the entire ^He 
phase space in the pd^^Heij reaction at these energies. The 
maximum emission angle of the ^He in the laboratory system 
is 18.5° at Tp=l36Q MeV and 19.6° at rp=1450 MeV and the 
FD only covers angles up to 18.0°. 

In figure [T] the acceptances at both energies are shown as a 
function of cosO*, when constraints optimised for rj ^ yy se- 
lection (see section [J. 1.1 1 ) are applied. The acceptance drops at 
high and low angles due to ^He ions emitted at small laboratory 
angles, ft/^^ < 3°. The middle hole in the acceptance is caused 
by^He ions emitted at large angles ft^^ > 18°. The acceptance 
drops at cos0* « -0.75 (1360 MeV) and cos6»* « -0.55 (1450 
MeV) are caused by ^He ions stopping between two layers of 
the FRH. 
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Table 1. Tiie constraints applied for selection of pd — > Hez/, 77 — > 77. 
The angle d(yyjn,m(^He) < 20° refers difference between the direction 
of the 77 system and the missing momentum of the ■'He. 

^He giving signal in the FPC and stopping in the FRH 
^ 2 photons in the SEC with Ey > 20 MeV 
one 77-combination fulfilling I /M(77) —m^\ < 150 MeV/c^ 
MM^i^Hejj) < 10000 {MeV/c^f 

^(yy)mm{^He) < 20° 

no overlapping hits in the FSB and the SEC 

160°<|fe,,(3He)-fc,,(77)|<200° 



3.1 Event selection 



The three main decay channels of the 77, i.e. rj-^yy (BR=39.3%), 

;r'';r'';r'' (BR=32.6%) and 77 ^ 7fn+n- (BR=22.7%) have 
all been separated and studied with the WASA setup. In this 
work, we focus on // — > yy since it provides a clean sample 
with good statistics. The simultaneous study of 77 7r'':7r'';r*' 
and 77 — > n^n^n^ allow valuable cross checks of the results. 



3.1.1 pd — > He77,77 yy 

In this case all final state particles - one ^He and two photons - 
can be measured with good acceptance. We thus have an over- 
constrained measurement and thereby, we can check if an event 
is consistent with the expected kinematics. This reduces the 
background significantly and gives a clean sample. 

The criteria for rj yy selection are given in table [1] As- 
suming phase space production, they give an acceptance of 
20% at 1360 MeV and 14% at 1450 MeV. 
Figure [1] shows how the acceptance varies as a function of 
COS0*. The acceptance is limited by the geometrical coverage 
of the FD, by photons missing the Csl modules in the calorime- 
ter and by the efficiency reduction due to ^He ions undergoing 
nuclear interaction before depositing all their energy. 

The upper panel of figure|2]shows the ^He missing mass for 
all events fulfilling the constraints optimised for t] ^yy selec- 
tion at rp=1360 MeV. The bottom panel shows the 7^=1450 
MeV case. Phase space Monte Carlo simulations of the main 
background channel, pd — > ^HcTr'^Tr'', are also shown, normalised 
to fit the data. They reproduce the background in the experi- 
mental data fairly well, except for an enhancement at high ^He 
missing mass at 7^=1450 MeV which is caused by 
pd ^Hew, a> n^y events that accidentally satisfy the cri- 
teria. Assuming phase space 2n'^ production give an acceptance 
of 3.6% at rp=1360MeV and4.0% at rp=1450MeV. Other re- 
actions, e.g. pd ^He;r*';r'':7r'', were found to give a negligible 
contribution to the rj ^yy background. 



3.1.2 pd ^ ^Hej],!]^ 7r%^7T° 

In this case we need six photons from the three n'^ decays in 
order to indentify the events. The Scintillator Electromagnetic 
Calorimeter (SEC) has a small "hole" in the backward part and 
one large in the forward part, where the photons escape un- 
detected. Therefore, in most 77 n^^jf'n^^ events at least one. 



so 



Acceptance at T^=1360 MeV 
Acceptance at T^=1450 MeV 
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Fig. 1. The acceptance of the pd — ^•'He77, J] ^ jy reaction as a func- 
tion of COS0* at rp=1450 MeV (dashed fine) and at 7^=1360 MeV 
(solid line). 
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Fig. 2. The upper panel shows the WASA data sample fulfilling the 
constraints optimised for selection of 77 —s- 77 at 1360 MeV and the 
lower panel shows the 1450 MeV case. The solid line histograms show 
Monte Carlo simulated pd "'He;7r"7r" data fulfilling the given con- 
straints. The spectra are not corrected for acceptance and the back- 
ground simulations are scaled to fit the data. 



but often several, photons escape detection. The acceptance is 
therefore significantly reduced compared to the rj ^yy case. 

The constraints optimised for 77 — s> selection are given 

in table |2] 

Assuming phase space production, this gives a total acceptance 
of 5.7% at rp=1360 MeV and 3.6% 7^=1450 MeV. The main 
background channel is direct pd -^^Weifiifiyfi production. At 
high missing masses, there is also a contribution from 

production, which will be discussed in sec- 
tion l4.1l The acceptance for direct 37r'' production at 7'p=1360 
MeV is 1 1.7% and 10.3% at rp=1450 MeV, if phase space pro- 
duction is assumed. 

The upper panel of figure[3]shows the ^He missing mass for 
all events fulfilling the constraints optimised for 77 n^n^yfi at 
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mesons in pd — s-^'HeX 



Table 2. Tiie constraints applied for selection of 
pd -^^ HcTy, 77 — ^ ;r";r";r". 



He giving signal in the FPC and stopping in the FRH 
> 6 photons in the SEC with Ey > 20 MeV 
one 77-combination fulfilling \IM{yy) — m^o \ < 50 MeV/c^ 
two other 77-combinations fulfilling |/M(77) — ra^o | < 60 MeV/c^ 
MM^{^lle6y) < 20000 (MeY/c^f 
no overlapping hits in the FSB and the SEC 




Fig. 3. The upper panel shows the data fulfilling the constraints opti- 
mised for selection of j] ^ rfiifiTfi at 1360 MeV and the lower panel 
shows the 1450 MeV case. The dotted line histograms show Monte 
Carlo simulated pd — s- ■'He;r data fulfilling the given constraints, 
the dashed-dotted histogram simulated pd — s- "'He;^ data and 

the solid line the sum of 3;^^ and Arfi production. The spectra are not 
corrected for acceptance and the background simulations are scaled to 
fit the data. 



Table 3. The constraints applied for selection of 
pd -^^ He?;, 77 — i> n'^n^n^K 

^*He giving signal in the FPC and stopping in the FRH 

> 2 photons in the SEC with Ey > 20 MeV 

one 77-combination fulfilling I /M(77) — nvl < 45 MeV/c^ 
MMi^n^yfi) > 250 MeV/c^ 

> 2 hits in the FSB 

£tot(SEC) <900 MeV 



rp=1360 MeV and the bottom panel shows the same but for 
rp= 1450 MeV. 




450 500 550 600 650 700 750 800 
MM3He(MeV/c^) 



Fig. 4. The upper panel shows all data at 7^=1360 MeV that satisfy the 
criteria optimised for pd —5- ^Hej],}] — s- n^n^jfi selection. The solid 
line represents Monte Carlo simulations of direct Tr+^r^^r" production. 
These spectra are not corrected for acceptance and the background 
simulations are scaled to fit the data. The lower panel shows the same 
thing but in the angular region 0.6< cos^?* <0.8. The line is the result 
of a fit of a gaussian peak on top of a polynomial background. 



The main background comes from nonresonant jt^tt^jt 
production. The acceptance for the pd — s> ^He7r+ tt^tt'' reac- 
tion when the given constraints are applied and phase space 
production is assumed, is 35% at Tp = 1360 MeV and 31% at 
Tp = 1450 MeV. 

The upper panel of figure|4]shows all data at Tp= 1 360 MeV 
that fulfill the cuts optimised for pd ~¥^Heri,ri ^ n^n^n^ se- 
lection. It is difficult to separate the rj events from the back- 
ground, partly due to the small signal-to-background ratio and 
partly due to the broad t] peak. However, in individual regions 
in cos6*, the 77 events appear in a peak and can be separated 
from the background with reasonable accuracy. An example is 
shown in the lower panel of figure |4] The q peak for the full 
cos9* range, shown in the upper panel of figure IH is broader 
than the 77 peak in an individual cos 6* interval, shown in the 
lower panel of figure |4] The broadness of the peak in the full 
cosO* range is due to a small dependence of the 77 peak position 
on cos 6*. This in turn is an effect of the calibration constants, 
which are slightly dependent on energy. This was also observed 
in Ref. Q, but there the effect was much stronger. Here it is 
negligible for small lab angles 0''^^ where the variation in Ts^^ 

is small. For large 0''^^, it gives a contribution to the systematic 
uncertainty of < 3 %. 



3.1.3 pd ^^Her],!] ^ TT^TT n'^ 

The criteria optimised for 77 — > 7r+ tt^tt" are given in table [3] 
The last one, requiring the total energy deposit in the SEC to 
be smaller than 900 MeV, rejects time-overlapping events, i.e. 
chance coincidences. The selection criteria give altogether an 
acceptance of 18% at Tp = 1360 MeV and 12% at Tp = 1450 
MeV. 



3.2 The 77 angular distribution 

The angular distributions were obtained by dividing the rj^yy 
data sample into intervals of cos0* where the acceptance is 
smooth and non-zero. The 77 mesons are identified by the miss- 
ing mass method in individual bins of cos0*. The 77 mesons are 
easier to identify in the intervals than in the cumulative spec- 
trum (compare the upper and the lower panel of figure |4]l. The 
number of 77 candidates is extracted by fitting Gaussian peak on 
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Fig. 5. The ■'He missing mass distribution for events satisfying tfie 
constraints optimised for rj ^ yy selection (see text) in given inter- 
vals of cos 9*. These MM-spectra are not acceptance corrected and 
the background simulations are scaled to fit the data. 



top of a polynomial background (it has been checked that in in- 
dividual cosO* region the background has no discontinuities). 
This number was then corrected for acceptance. The system- 
atic uncertainty was estimated by fitting simulated Monte Carlo 
data of the main background channel (in this case 
and compare the number of 77 events obtained in this way to the 
number of 77s obtained from fitting the background to a polyno- 
mial. The same procedure was repeated for the 77 3 tt chan- 
nels. It turns out that the agreement in individual cos 6* regions 
is good between the t] 2j and the 77 —s> 3 tt channels. This 
gives confidence that the cut efficiencies are well understood 
and that our systematic uncertainties are under control. 

The normalisation was achieved by comparing data on back- 
ward going 77 mesons from pd ->3He77 from SPES IV (U) and 
SPES III [12] with the corresponding data from this work using 
the method described in Refs. 1231 and Il25i l. The normalisation 
uncertainty of the measured cross sections and is 29% at 1360 
MeV and 12% at 1450 MeV. 

The resulting angular distributions are shown in figure |6] 
and figure|7] The systematic uncertainties are shown as a shaded 
histogram in each figure. They mainly arise from the ambigu- 
ity in the background subtraction, but there is also a small con- 
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Fig. 6. Angular distribution of the 77 meson in the CM system at 
Tp = 1360 MeV. The black dots are WASA data from the pd ^■ 
"^He;;,;; —5- 77 channel. The error bars represent the statistical uncer- 
tainties and the grey histogram the systematical. The open dots are 
obtained with WASA data from pd — > 'He77,77 — ?■ ;r'';r+;r^. The black 
triangles are calculated by interpolating SPES IV 1 1 1] data at r^=1250 
MeV and 7^=1350 MeV and SPES III data at 7^=1450 MeV. The 
Legendre fit represented by the solid line has been used to calculate 
the cross section quoted here and fit shown as a dashed line has been 
used to estimate the systematical uncertainty. 



tribution from the energy dependence of the calibration con- 
stants (see section 13.1.3b . The distributions at both energies 
are highly anisotropic with a sharp forward-backward asym- 
metry. This is in Une with earlier experiments, e.g. Refs. ||7][8] 
l9l[T0l[T2l. where evidence were found for several higher partial 
waves away from the threshold region. From comparing SPES 
III data with data from this work at rp=1450 MeV, which is 
done in figure |7] the conclusion is that either the two data sets 
are inconsistent, or there is a forward dip that is much stronger 
than the dip observed in Refs. ||7][8l. 

The angular distributions were fitted by a series of Legen- 
dre polynomials 



52 akPk{cose*) 

k=0 



(4) 



to the T] ^ yy data points from WASA. The zeroth coefficient 
of the Legendre polynomial gives, when multiplied with 47r, 
the total cross section. At 1360 MeV one obtains 
o"tot = 151.6 ±9.3 ± 35.3 nb. In addition, there is an uncertainty 
from the normalisation of 29%. At 1450 MeV the total cross 
section is estimated to be ixtot = 80.9 ± 3.6 ±43.0 nb. The nor- 
malisation uncertainty is 12% at 1450 MeV. 



4 Multipion production 

In this section, we first study the pd — ^^He7r'';r''7r'' reaction, 
then the pd ^>^He7r+ tt^tt" reaction and finally, the two three- 
pion reactions are compared. 
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Fig. 7. Angular distribution of the rj meson in the CM system at Tp = 
1450 MeV. The black dots are WASA data from pd •'He?/, rj^yy 
at 7^=1450 MeV. The error bars represent the statistical uncertainties 
and the grey histogram the systematical. The open dots are obtained 
with WASA data from pd ^ller],r] Tfin^n^ . The black triangles 
are data from SPES III 1 12| while the open triangles come from SPES 
rV nil . The curves are results of fits of Legendre series to the WASA 
data (solid) and WASA plus SPES III (dashed). 

4.1 The pd ^^HejT^TT^TT^ reaction 

The same selection criteria are used as for the pd ^-^He ri,ri ^ 
ifin^jfi case, given in table |2] of section U. 1.21 For the 
pd ^-^HciTrVTr" reaction, they give acceptances of 11.7% at 
1360 MeV and 10.3% at 1450 MeV. There may be a large un- 
certainty in the acceptance of a reaction where six photons are 
measured. To estimate this uncertainty, we assume that the dif- 
ference in the extracted number of 77 mesons from rj ^yy and 
7] — !■ 7r'^;7r'';7r'' is entirely caused by the ambiguities in the accep- 
tance and that the uncertainty is the same at both energies. The 
uncertainty in the acceptance is then estimated to a maximum 
value of 20%. 

The ^He missing mass distributions at both energies for all 
events fulfilling the constraints are shown in figure[3]in section 
13.1.21 The dotted line shows simulated pd ^HeTr^Tr^Tr'^ data 
assuming phase space production. Simulated 37r'' data match 
the experimental data for low and medium missing masses (ex- 
cept at the r] peak, which is expected), but at high MM(^He), the 
matching between data and phase space Monte Carlo is poor. 

It is reasonable to assume a contribution from 

either from direct production or from pro- 
duction via t/tt" in pd ^HeT/Tr",;; — s> Tr'^Tr'V^. In both reac- 
tions, eight photons are produced and the acceptance for the 
selection criteria in table |2] is 28% at 1360 MeV and 24% at 
1450 MeV. At the highest energy, the maximum ^He emission 
angle in the lab system is 15° in the rin^ case and 18° in the 47r" 
case, which means that in both cases, the WASA Forward De- 
tector covers almost the full ^He phase space. The acceptance 
is then nearly independent of the production mechanism. 

The 4-7t'^ distributions obtained from Monte Carlo simula- 
tions are shown in the dashed-dotted line histograms in the up- 
per and lower panel of figure |3] Adding the contributions from 
37r'' and A-n^ together gives the solid line histograms in figure 



|3] We obtain N^^o = 250 at Tp = 1360 MeV and N^^o = 800 at 
Tp = 1450 MeV. 

The cross section of the pd — > ^Hez/Tr*^ reaction at Tp = 
1450 MeV has been measured to (Ttot = 23.6 ± 1.6 ± 2.2 nb 
±14% by studying the rj ^yy decay lIZTl . 

From the known cross section of pd ^He 77:77'' reaction 
at Tp = 1450 MeV (see Ref. EtI ) and from the acceptance 
and the branching ratio of 77 ^> 7r'';7r''7r'', the expected number 
of pd ^He.T]7T^,i] n^TJpjfi events is calculated to 700 ± 80. 
This explains almost fully the A^4^o = 800 and it is clear that 
the cross section of direct 4;7r'' production must be very small at 
1450 MeV. 

Subtracting the fitted AtiP and rjifi distributions from the 
experimental data gives A^3^o = 1400 and A^3^o = 4500 at Tp = 
1360 MeV and Tp = 1450 MeV, respectively. This coiTesponds 
to 3?:'' cross sections of 180 nb and 115 nb. The statistical un- 
certainty of A^3jj.o is given by the square root of the total number 
of events before the subtraction and is equal to 3% (2%) at 
rp=1360 MeV (rp=1450 MeV). 

It is also possible that part of the deviation from the 377^ 
phase space curve in figure[3]is due to a production mechanism 
that differs from phase space production. This will be discussed 
later in this paper. However, at least at rp= 1450 MeV, the ex- 
pected contribution from pd — ?► ^HcT/jr", i] 2>n^ explain the 
data well and the remaining excess of events at high ^He miss- 
ing masses gives a small contribution to the systematic uncer- 
tainty. 

At rp=1360 MeV, it is difficult to say with certainty that the 
excess of events at high MM(^He) in the upper panel of figure 
[3]are not directly produced ^jfi events. The cross section of the 
pd^^We. rjn^^ reaction is not known, and it is therefore unclear 
whether a significant contribution from this reaction is to be ex- 
pected. However, the mixture of and An^ events reproduces 
the experimental distributions also at Tp = 1360 MeV very well 
and it is therefore reasonable to assume a contribution from 
471" production, either from direct production or from the sub- 
sequent 77 decay in rijf' production. We therefore take the 3?:'' 
cross section of 180 nb, calculated when assuming that the de- 
viation from the ^jfi curve at large MM(^He) in figure[3]come 
from 4;7r*' production, as the most reliable one. The excess of 
events is treated as a systematic uncertainty. By assuming that 
all events in the upper panel of figure [3] that do not come from 
pd — >^He rj,?]^ 37:'^ are directly produced 37:'^ events, N^^o be- 
comes 1650 which coiTesponds to a cross section of 212 nb. 
The systematic uncertainty is then taken as the difference be- 
tween the cross sections calculated in two different ways, i.e. 
32 nb. We assume that the uncertainty is symmetric. This is a 
conservative method of estimating the systemtaic uncertainty 
and other systematic contributions, e.g variation in the accep- 
tance due to reaction mechanism, should be well within the 
error bars estimated in this way. 

We can also give a rough upper limit of the pd — >^'He77;7r'' at 
1360 MeV, which will be useful in the next section. Assuming 
that all the N^j^o = 250 events come from 7771" production, the 
o-{pd -^^HcT/jr") would be 42 nb. 

The total cross section of 37!-'* production then becomes 
cTg^o = 180±6±49 nb ±29% at Tp = 1360 MeV and 0-3^0 = 
1 15 ± 3 ± 23 nb ± 12% at = 1450 MeV. The first error is sta- 
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Fig. 8. The missing mass of the ■'He;r"-system, which is equivalent to 
the invariant mass of the 27r''-system. The points represent background 
subtracted and acceptance corrected data satisfying the criteria given 
in the text and 600 MeV/c^ < MM(^He) < 700 MeV/c^. The sohd 
histogram shows phase space Monte Carlo simulations of StjP produc- 
tion. The upper panel shows the Tp = 1360 MeV case and the lower 
the Tp = 1450 MeV case. 




3500 



IM.„^,. (MeV/c^) 

Fig. 9. The invariant mass of the ■'He7r*'-system. The points repre- 
sent background subtracted and acceptance corrected data, fulfilling 
the constraints given in the text and 600 MeV/c^ < MM(^He) < 700 
MeV/c^. The solid histogram shows phase space Monte Carlo simu- 
lations of Tiifi production. The upper panel shows the Tp = 1360 MeV 
case and the lower the Tp = 1450 MeV case. 



tistical, the second is systematical and includes uncertainties 
from background and acceptance. The last uncertainty comes 
from the normalisation. 

Background from quasi-free reactions pp — > ppn'^n^n'^ with 
a proton misidentified as a ^He is expected to be negligible. The 
probability that an event from a reaction with p or d in the final 
state instead of ^He would survive the constraints, is smaller 
than 0.001%. 

Invariant mass distributions of the final state particles in the 
pd — ^^HeTT^TT^Tr" reaction give important information about the 
production mechanism. Deviation from phase space can give 
hints about e.g. intermediate resonances. In this work we have 
studied the 2:7r''-system, the ^HeiTr^-system and the ^He2;7r''- 
system. When studying the invariant mass of two pions it is 
more convenient to instead reconstruct the missing mass of the 
^He and the third pion, here denoted MM(^Hen'^). This is be- 
cause the ^He is measured in the FD with higher resolution than 
the pions, which are measured in the CD. 

In order to avoid an event sample with a lot of background 
from the pd -^^Hej] and pd -^^Herin'^, events which fulfill 
the condition 600 MeV/c^ < MM{^He) < 700 MeV/c^ are se- 
lected. MM{^Hen'^) is then reconstructed for these events. In 
this event sample, there will be a small contribution (a few per- 
cent) from the pd ^^He ijtt^, rj — > Stt" reaction and the data will 
therefore be subtracted by the expected amount of rjn'^, rj irfi 
events, obtained from simulations. The data are then corrected 
for acceptance. The results are shown in the upper and lower 
panel of figure [8] The points represent the background sub- 
tracted and acceptance coiTected data and the solid histogram 
phase space simulated 2>n^ data. The experimental data follow 
phase space well. 

The invariant mass of the ^He:7r''-system, /M(^He:7r''), has 
also been reconstructed. The small background from rjn^ was 
subtracted in the same way as in the MM(3He7rO) case and the 
data was then corrected for acceptance. The result is shown in 



figure |9] Here the data disagree with phase space. There is an 
enhancement with respect to phase space centered around « 
3090 MeV/c^, which roughly equals the sum 2mp + M^^^y^ii)- 
This may indicate a single zl( 1232) excitation in the production 
process. 

Finally, the invariant mass of the ^He7r'':7r''-system is stud- 
ied. This is done by reconstructing the missing mass of one of 
the tt" mesons. In this way, the resolution is improved. The ex- 
perimental data were background subtracted, using simulated 
pd -^^Werju^ ,7] data, and acceptance corrected. The re- 

sult is shown in figure[TOl together with the phase space Monte 
Carlo simulations of pd — s-^Hctt There is a small en- 

hancement at high MM{7f') around the Inip + M^«(i44o) sum, 
which may indicate the involvement of a Roper A^*(1440) ex- 
citation in the production mechanism. 

4.2 pd^'^ He;r+;r-;r'' 

All selection criteria optimised for 77 — > tt+tt^tt" selection, given 
in table [3] in section !?. 1.31 are applied. In addition, we require 
that both charged pions are emitted in directions covered by the 
CD, i.e. that no other charged tracks than the ^He are found in 
the FD. Finally, two charged tracks in the MDC are required, 
with overlapping hits in the PSB. The total acceptance, cal- 
culated partly using Monte Carlo simulations (all constraints 
not involving the MDC, see section [J. 1.31 ) and partly analysing 
cj — !■ -K^n^n^ data (all constraints involving the MDC, see 
Ref.|23|) is 7.2% at 1360 MeV and 6.7% at 1450 MeV. 

In figure [TT] the missing mass of the ^He is shown for aU 
events satisfying the criteria given in table |4] There is a small 
enhancement around the r] mass and a clear peak at the w mass, 
but except from that, the experimental data seem to follow the 
phase space ;7r+;7r^7r" distribution well. There is no sign of any 
pd -^^W&rpi^ ,1] n^n^n^ events in this sample. The accep- 
tance for this reaction when applying the cuts in table |4] is 9% 
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Fig. 10. The missing mass of the jfi 
ant mass of the ' 

tracted and acceptance corrected data fulfilling the constraints given 
in the text and 600 MeV/c^ < MM(^He) < 700 MeV/c^. The solid 
histogram shows phase space Monte Carlo simulations of S/r" produc- 
tion. The upper panel shows the Tp = 1360 MeV case and the lower 
the Tp = 1450 MeV case. 
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Fig. 11. The upper panel shows the WASA data sample fulfilling the 
constraints optimised for selection of pd — s-^ Yle-n^n^rfi at 1360 MeV 
and the lower panel shows the 1450 MeV case. The solid line his- 
tograms show Monte Carlo simulated pd — > ■'He7r+;r^;r*' data fulfill- 
ing the given constraints. The peak in the experimental data at high 
missing masses are pd — s--' Heoj, w — i- Tt^n^iP events. The spectra 
are not corrected for acceptance and the background simulations are 
scaled to fit the data. 



Table 4. The constraints applied for selection of pd — s--' He7r+7r 

^He giving signal in the FPC and stopping in the FRH 
^ 2 photons in the SEC 

one 77-combination fulfilling \lM{yy) — rrijf, \ < 45 MeV/c^ 
MM(3He;r°) > 250 MeV/c^ 
^ 2 hits in the FSB 
£tot(SEC) < 900 MeV 
no in FD 

2 MDC tracks with matching hits in the PSB 



at 1360 MeV and 12% at 1450 MeV. At the higher energy, 
where the t^tt'^ cross section is known (see [[271 ). the expected 
number of pd -^^W&rpi^ ,t] — > n^n^n^ events is «i250, which 
constitute only w 1 % of the continuum data in the lower panel 
of figure [TT] In the previous section we found that the direct 
AtP cross section at 1450 MeV must be very small, and even 
if the cross section of direct n^Ti^Tx^Tfi is likely higher than 
the direct 4;r*' cross section due to more possible isospin am- 
plitudes, it is reasonable to assume that direct n^n^Tf'Tf' pro- 
duction will give a negligible contribution to our n^n^if' data 
sample. This assumption is also very well in line with the good 
agreement between data and simulations in figure[TT] The num- 
ber of n^n^if' events at 1360 MeV, A^3;r = 6700, corresponds 
to a total cross section of 1400 nb. This is obtained using an 
acceptance which is calculated assuming phase space produc- 
tion, but since the WASA detector covers the major part of the 
^He phase space for this reaction, the model dependence of the 
acceptance is small. 

The largest contribution to the systematic uncertainty comes 
from the efficiency of the MDC. A robust method of estimating 
this uncertainty is to calculate the cross section with and with- 
out using the information from the MDC, treat the difference 
as a systematical uncertainty and assume that it is symmetric. 



At 1 360 MeV, the number of ii^n^n events obtained using 
selection criteria without involving the MDC is N^„ = 33000 
and the acceptance is 28%. This corresponds to a cross section 
of 1770 nb which gives a systematical uncertainty of 370 nb. 

There may be a systematical uncertainty arising from falsely 
identified pd —^^¥lej]K^,ri — > 7r+ tt^tt". According to the rough 
upper limit of the t/itt" cross section at 1360 MeV that was given 
in the previous section, the maximum number of t/tt" events in 
this sample is 56. This gives a systematic uncertainty of 12 nb 
and is thus very small compared to the uncertainty from the 
MDC efficiency. 

The statistical uncertainty is obtained by the square root of 
the total number of events and is determined to 17 nb. Finally 
we get o-„^„-„o = 1400 ± 17 ±370 nb ±29% at Tp = 1360 MeV. 

Following the same reasoning at 1450 MeV, except that the 
rjn^ cross section is known with good precision, the cross sec- 
tion becomes cr^+^-^o = 910±7±80nb ±12%. 

As in the 3n^ case, the distributions of the final state par- 
ticles have been studied. Since the tt" is the only pion that 
is fully reconstructed, the invariant mass of the 7r+;r^-system 
is studied by reconstructing the missing mass of the ^HeTr*'- 
system. To have a sample as similar to the 3;r'' case as possible, 
events which satisfy 600 MeV/c^ < MMC'Hc) < 700 MeV/c^ 
are selected. The MM{^HeK^) is reconstructed and the data 
are corrected for acceptance. Note that no background subtrac- 
tion had to be made in this case, since the contribution from 
pd — !>^He777r" is proven to be small at both energies. The re- 
sults at both energies are shown in figure [12] There is good 
agreement between experiment and simulated n^n^n^ data and 
there is no sign of any intermediate p meson, which would push 
the MM{^Hen^) towards higher masses. This is not surprising 
since despite the large width of the p meson (/" « 150 MeV), 
we are far below the nominal pd ^^^Hep'^Tr" threshold at both 
beam energies considered in this work. 
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Fig. 12. Thie missing mass of the •'He;n"^-system. The points represent 
acceptance corrected data satisfying the criteria given in the text and 
600 MeV/c-2 < MM{^H&) < 700MeV /c^. The solid histogram shows 
phase space Monte Carlo simulations of n'^n^ifi production. The up- 
per panel shows the Tp = 1360 MeV case and the lower the Tp = 1450 
MeV case. 
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Fig. 14. The missing mass of the ;r". The points represent accep- 
tance corrected data fulfilling the constraints given in the text and 
600 MeV/c2 < MM(^He) < 700 MeV/f^. The solid histogram shows 
phase space Monte Carlo simulations of 7r+7r^;r" production. The up- 
per panel shows the Tp = 1360 MeV case and the lower the Tp = 1450 
MeV case. 




(MeV/c^) 

Fig. 13. The invariant mass of the ^He/r'^-system. The points represent 
acceptance corrected data fulfilling the constraints given in the text 
and 600 MeV/c^ < MM(^He) < 700 MeV/c^. The solid histogram 
shows phase space Monte Carlo simulations of n^jr^n^ production. 



The upper panel shows the Tp ■ 
= 1450 MeV case. 



1360 MeV case and the lower the T,, 



The invariant mass of the ^HeTr" is reconstructed in the 
same way as in the Stt" case, except that no background sub- 
traction was necessary. The results are shown in figure[T3] Also 
in the n^n^n^ case, there is a small enhancement around the 

2m„ + M, 



zl(1232) 



sum. 



The invariant mass of the ^He.7r+7r system is studied by 
reconstructing the missing mass of the n'^. The acceptance cor- 



rected data are shown in figure[T4] Like in the pd -^^He7r'V^;7r' 
case, there is a small enhancement with respect to phase space 
near the Inip + M^«(i44o) sum, suggesting that the Roper reso- 
nance A^* (1440) may be involved in the production process. 



4.3 Comparison between the pd 

pd — He7r''7r''7r'' reactions 



HCTT^ 



tt" and the 



In the introduction, the ratio between the cross sections of the 

pd —5 

(r{pd- 



HeTT+TT tt" and the pd -^^ HeTr^Tr^Tr'' reactions, i.e. 



, , -vT, — n-rnrr, was calculated to 9, using a statistical model 
where all isospin amplitudes Mj^^ are put equal and all cross 
terms are set to zero. If instead Mq = 0, the ratio becomes 4. In 
this work, the ratio has been measured experimentally at both 
energies. By using the cross sections determined in section 143] 
andgJl one then obtains 7.8 at 1360 MeV and 7.9 at 1450 
MeV for this ratio. However, to give a comparison at the same 
excess energy Q, the results have to be corrected for the differ- 
ence between the masses of the tt* and jr". The lower mass of 

.o„o„o 



^Hctt' 

^ He 7r+ tt^tt'^ at the same 



the tt" makes the phase space volume of the pd 
reaction larger than that of the pd 
beam energy. After correcting for the difference in phase space 
volume, the ratio becomes 8.3 ±0.3 ±3.1 at 1360 MeV and 
8.4±0.2 ± 1 .8 at 1450 MeV, where the first uncertainty is statis- 
tic and the second systematic. Note that the uncertainty in the 
normalisation cancels in the ratio. The values obtained are con- 
sistent with the value of 9 predicted using the statistical ap- 
proach. The interpretation of this result is then that Mq should 
be of similar size as Mi . 



5 Summary and Conclusions 



The production of light mesons, i.e. n and 77, have been stud- 
ied at Tp = 1360 MeV and Tp = 1450 MeV. The pd ^^He/? 
reaction was studied by using data from the three most com- 
mon decay channels; i] yy, rj n'^n'^n'^ and rj — > :7r+;7r^;r'^. 
The result from the different channels gave consistent results. 
At both energies, the angular distributions of the 77 meson were 
reconstructed and they show a pronounced forward-backward 
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Table 5. The total cross sections of the reactions studied in this work. 
In addition to the systematic uncertainty given in the table, there is 
a normalisation uncertainty of 29% at 1360 MeV and 12% at 1450 
MeV. The normalisation is made using pd -^^Hej] data and differen- 
tial cross sections given in I11I12I . 



Reaction 


Tp (MeV) 


cr 


Stat. 


syst. 






inb) 


inb) 


(nb) 


pd -^^Heri 


1360 


151.6 


±9.3 


±35.2 


pd ^3He;r°;r°;r° 


1360 


180 


±6 


±49 


pd -^^Hen+n-n^^ 


1360 


1400 


± 17 


±370 


pd -^^Hei] 


1450 


80.8 


± 3.6 


±43.1 


pd -^^Hen^n^n^ 


1450 


115 


± 3 


±23 


pd -^^Yien^n^TiP 


1450 


910 


± 7 


±80 



asymmetry. The WASA detector does not cover the full angu- 
lar range at these high energies and one could therefore not say 
whether the forward plateau or dip observed in |7,8| persists 
at high energies. The data from this work in combination with 
data taken at the same energy with the SPES III spectrometer 
lfT2l do, however, suggest that the forward dip persists at this 
high energy and is stronger compared to the lower energy case. 
If not, the 1450 MeV data from this work disagree with the 
SPES III data in the forward hemisphere. 

The total cross sections of three pion production in 

and pd ^-^HcTr+Tr tt" were measured at both 
energies. The ratio between the n'^n^n'^ and Stt" cross sections 
was calculated at both energies and the results are consistent 
with the statistical model, where Mq — Mi and all cross terms 
are neglected. The invariant mass distributions of the two-pion 
system and the ^Hctt" system were reconstructed. The invariant 
mass distributions of the two-pion system follow phase space 
but the corresponding distribution of the ^HeiTr" system, show a 
small enhancement around the 2mp + ^^(1232) mass. This en- 
hancement was observed in 3n'^ and n'^ji^n'^ production at both 
energies and may indicate a single A excitation in the produc- 
tion mechanism. 

The invariant mass distributions of the ^He2;7r system, for 
both n^n^n'^ and 3n'^ production at both eneregies, show an 
enhancement near the Imp + M^*(i44o) mass, suggesting that 
the Roper A^* (1440) resonance may be involved in the produc- 
tion mechanism. 

The cross sections measured in this work are summarised 
in table 21 
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